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1 Introduction 

The spectrum of baryon resonances is expected to be very 
rich and considerably more complex than the mesonic 
excitation spectrum. Yet, experimentally, the number of 
known light-quark mesons exceeds by far the number of 
known baryon resonances [1]. 

In quark models [2-4], most high- mass baryon reso- 
nances are only weakly coupled to Nir [5], and can thus 
not been seen in elastic ttN scattering experiments. For in- 
elastic reactions like irN -> rjN, irN -> KA, irN — > K£, 
data with a polarized target are missing, and data on 
differential cross sections have low statistics, and are of- 
ten inconsistent when different experiments are compared. 
States weakly coupled to the ttN channel may thus have 
escaped identification. The situation was aggravated by a 
recent analysis of a large body of irN elastic and charge 
exchange scattering data in which many of the less estab- 
lished nucleon and A resonances were not confirmed [6]. 

Other interpretations of the baryon spectrum exist as 
well. Very popular are diquark models [7-10] in which one 
quark-pair is frozen and in which the number of predicted 
states decreases. Further, we mention approaches based on 
chiral Lagrangians in which low-lying baryon resonances 
are generated dynamically. In many cases, these calcula- 
tions offer a consistent description of resonance properties 
and scattering data (see e.g. [11]) but so far, they do not 
give a survey of all resonances to be expected. Often dis- 
cussed is the conjecture that chiral symmetry might be 
restored in high-mass meson and baryon resonances. [12, 
13]. The conjecture gives an attractive interpretation of 
one experimental observation, that resonances show up as 
parity doublets or even higher multiplets. It fails to pre- 
dict at which mass resonances should be found and, exper- 
imentally, all meson and baryon resonances on the leading 



Regge trajectory have no parity partner. The AdS/QCD 
model describes QCD in terms of a dual gravitational the- 
ory [14, 15]; with some phenomcnological adjustments, it is 
surprisingly successful in predicting the baryon mass spec- 
trum and the number of expected states [16,17]. It also 
predicts where parity multiplets should occur and where 
not. Reviews of baryon spectroscopy can be found in [18- 
20]. 

A decision which of the above approaches provides the 
most accurate representation of Nature requires a bet- 
ter experimental knowledge of the excitation spectrum. 
The limitations of the current data base on light-quark 
baryons is the stimulus for experiments studying baryon 
resonances in photoproduction of complex final states where 
the irN channel can be avoided in both the initial and 
the final state. The analysis of multibody final states in- 
cluding fermions is complex; in order to compare different 
approaches and to identify possible problems, a common 
meeting ground is needed. This is provided by the simplest 
photoproduction reactions, by jp — > pir° and jp — > nir + . 
In this paper, we give masses, widths, photocouplings, and 
Ntt decay branching ratios for the most important con- 
tributing resonances and compare our pion photoproduc- 
tion and helicity amplitudes to those obtained by SAID 
[21], MAID [22] and within the Giefien model [23,24]. The 
fits are based on a large number of data sets and include 
data with multibody final states. The study thus shows to 
which extent multibody final states are compatible with 
the best-studied Nir system. 



2 Data used in the fits 

A large number of reactions is used in the truly coupled- 
channel fits presented here. The data cover elastic ttN 
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Table 1. Pion induced reactions fitted in the coupled-channel 
analysis and x contributions. 
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scattering as well as inelastic reactions, they cover differ- 
ential cross sections and single and double polarization 
variables. Reactions with multi-body final states are in- 
cluded exploiting an event-based likelihood method. 

Different data sets often have a very different statis- 
tical power. Weights Wi are introduced to force the fit to 
take into account highly significant but low-statistics data, 
e.g. beam asymmetries. Without these weights, polariza- 
tion data often have too small an impact on the fit result. 
The weight of a newly introduced data set is increased 
when the fit is visually unacceptable, or decreased until 
first discrepancies between data and fit become apparent. 

2.1 Elastic tvN — > tvN scattering 

In the analysis presented here, data on elastic irN scat- 
tering (charge exchange is implicitly included) are not 
used directly. Instead, we rely on the detailed work of the 
George- Washington Center for Nuclear Studies [6] and use, 
for energies up to 2.2 GeV, their scattering amplitudes. 

2.2 The reaction iv~p -qn 

The inelastic n~p scattering process leading to the mj fi- 
nal state was reported from several experiments [25-28] 
and [29,30]. Above 1.8 GeV, large discrepancies between 
the data [27,30] show up. The data from [29] cover ex- 
treme backward angles and are partly incompatible with 
all other results. A critical discussion of the available data 
can be found in [31]. We use here the data from [26,28] 
(see Table 1) which show better consistency. 

2.3 The reaction 7r~p — > n 7v°n 

In the low-energy region, up to ~ 1.5 GeV in mass, very 
precise data from BNL are available [32]. These data are 
included in an event-based likelihood fit (see Table 2). 



Table 2. Reactions leading to 3-body final states are included 
in event-based likelihood fits. The xV-^bm values are calcu- 
lated from selected Dalitz plots (see text for details). Refer- 
ences to the data are given in the text. 
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The likelihood values have no direct significance; only like- 
lihood difference can be related to probability changes 
when particular contributions are removed from the fit. To 
demonstrate the quality of the description we have con- 
structed for every energy of the initial pion rnz n0 versus 
tt? o Dalitz plots for data and for Monte Carlo events with 
40 x 40 bins. The Monte Carlo events were weighted with 
the squared amplitude from our final PWA solution. The 
data and weighted Monte Carlo Dalitz plots were com- 
pared; in Table 2 the x 2 l^bin is given as well as the num- 
ber of bins with nonzero number of Monte Carlo events. 



2.4 Photoproduction of single neutral pions off protons 

References to the data on the reaction 7^ — > pir° and 
their x 2 contributions are collected in Table 3. For the 
differential cross section, we use only the most recent data, 
reported by TAPS@MAMI [41], GDH-A2 [42,43], GRAAL 
[44], CB-ELSA [45,46], and CLAS [47] which cover a wide 
range of energies and angles. A large variety of older data 
exist which cover only a limited fraction of the energy and 
angular range. These data provide significant information 
on polarization observables. 

The differential cross sections reported by the different 
collaborations exhibit small but significant systematic dis- 
crepancies practically in all mass regions; due to the small 
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Table 3. Observables from 7r and r\ photoproduction fitted in 
the coupled-channel analysis and x 2 contributions. For pion 
production, free normalization factors and additional system- 
atic errors were introduced to allow for data variation beyond 
statistical expectations (see text). 
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statistical errors these are easily recognized. We show the 
systematic deviations by comparing the data with a curve 
representing the "first" fit to all data, without normaliza- 
tion factors. 

In Fig. 1 the differential cross section from GDH-A2, 
CLAS and CB-ELSA are shown and compared to the 
preliminary fit for the 1495-1530 MeV mass range. The 
GDH data systematically exceed CLAS data while the 
CB-ELSA data provide numbers between these two mea- 
surements. The GDH-A2 data have a larger statistical er- 
ror than the CLAS data; we introduced them into the 
fit with larger weight since they provide important infor- 
mation about the difference between helicity 3/2 and 1/2 
cross sections. 

In the mass region 1600-1750 MeV, there are notable 
discrepancies between GRAAL and CLAS data (here the 
CB-ELSA data fall again between GRAAL and CLAS re- 
sults). As an example, the mass region around 1670 MeV 
is shown for the three data sets in Fig. 2. The curve cor- 
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Fig. 1. The CB-ELSA, CLAS and GDH differential cross sec- 
tion on 7p — > 7r°p in the region 1500 MeV. In the Crystal 
Barrel data, a common systematic error due to uncertainties 
in the reconstruction efficiency is included. The curve repre- 
sents the "first" fit without normalization (see text). The GDH 
data are introduced with a large weight. 
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Fig. 2. Comparison of three data sets on the jp — > u p differ- 
ential cross section in the region 1670 MeV. The curve repre- 
sents the "first" fit without normalization (see text). 



responds to the fit where the CLAS data are taken with 
statistical errors only and dominate the solution. 

At higher energies, only the CB-ELSA and CLAS data 
are available. There are two clear discrepancies between 
these data sets. The first one is located in the 1900 MeV 
mass region where the CB-ELSA data systematically ex- 
ceed the CLAS data in the backward hemisphere (see 
Fig. 3, top). A second discrepancy shows up above W = 
2100 MeV in the very forward angular range. Here the 
corresponding CB-ELSA points are systematically lower 
than those from CLAS. At IV=2300-2400MeV, the two 
data sets are fully consistent (see Fig. 3, bottom). 

A fraction of the discrepancies is assigned to normal- 
ization. Most experiments give explicit normalization er- 
rors which were not taken into account in the "first" fit. 
We then allowed for a free normalization factor for the 7r° 
differential cross section, which is determined in the fit to 
1.00 (TAPS@MAMI), 1.01 (GDH-A2), 0.99 (GRAAL and 
CB-ELSA) and 0.95 (CLAS). 

The data from the four experiments are still not yet 
statistically compatible; at least one experiment must have 
additional unrecognized systematic errors. Of course, we 
do not know which experiment. We assume that all four 
experiments have systematic errors which were not rec- 
ognized. These were estimated from the variance of the 
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Fig. 3. The CB-ELSA and CLAS differential cross section on 
7P -)• tt°p in the region of 1900 (top) and 2350 (bottom) MeV. 
The curve represents the "first" fit without normalization (see 
text). 



experimental results in preset bins of energy and angle (in 
cos 9). For this purpose, differential cross sections were 
calculated, by interpolation, for these bins. From the vari- 
ance we estimated systematic errors which increase lin- 
early from 1% at W = 1400 MeV to 9% at W = 2450 MeV. 
These systematic errors were added to all four data sets. 
With the increased systematic errors, the data are com- 
patible and the x 2 of a fit reflects the quality of a fit and 
not the inconsistency between different data sets. These 
errors arc used in the fits only. In the figures, the data and 
their errors are shown as quoted in the original papers. 

The beam asymmetry £ has been determined in a 
number of experiments [44,48-55], as well as the target 
asymmetry T [49-51,56-64], and the polarization P of 
the recoiling proton [49-51,65-69]. Few data exist from 
experiments with polarized photons and polarized target 
or from measurements of the recoil polarization. Data on 
O x and O z can be found in [68,72], on G in [70,71], and 
on H in [70]. Data on the helicity difference a 3 / 2 — <J\ii 
were published in [42,43]; in Tables 2 and 3 we quote E 
which is defined as (o'3/2~cri/2)/( cr 3/2 + f 7 i/2)- These data 
are included in the fits. Their statistical errors are mostly 
large, the systematic errors likely small. Hence we retain 
the original errors. 
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Fig. 4. Older data from different experiments (top) and CLAS 
differential cross section on 7p — > ir + n in the region 1630 MeV. 
The curve represents the "first" fit without normalization (see 
text). The consistency is excellent. 



2.5 The reaction jp — >• 7r+n 

Total cross sections were reported in [73-87] . Again, some 
discrepancies show up, at energies above 1600 MeV and 
in the forward region, between the new CLAS data and 
former measurements. An example of such discrepancies is 
given in Fig. 4. Normalization factors for the different data 
were introduced which arc determined to be in the range 
from 0.96 to 1.03. A consistent description was achieved 
by adding a systematic error which increases linearly from 
1% at W = 1400 MeV and to 9% at W = 2450 MeV. 

The beam asymmetry £ was determined in [55,88-98], 
the target asymmetry T in [93,94,99-109], the neutron 
recoil polarization can be found in [93,94,110]. A few data 
from double polarization are available: on G [71,111,112], 
H [111-113], and on the helicity difference a 3 / 2 — <J\/2 
[43,86]. The data are fitted with the errors as given in the 
respective papers. 



2.6 Photoproduction of r/ mesons off protons 

For photoproduction of rj mesons, differential cross sec- 
tions [114-117,119] and the related beam asymmetry £ 
[118-120] are the only quantities which have been mea- 
sured so far. Double polarization obscrvablcs arc presently 
studied intensively at several laboratories but so far, no re- 
sults have been published. The recent high-statistics mea- 
surements on 7p — > pr] [121,122] are not yet included in 
the fits presented here. 



2.7 The reactions -yp -> K+A,K+£° and K°S+ 

Data on hyperon photoproduction used in the present fits 
are collected in Table 4. We use the differential cross sec- 
tions for 7p -> K+A and K+£° from CLAS [123]. As 
shown in [124], the Saphir data [125] on differential cross 
sections are about compatible with the CLAS data when 
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Table 4. Hyperon photoproduction observables fitted in the 
coupled-channel analysis and x 2 contributions. 
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an energy dependent normalization factor is introduced. 
The beam asymmetry was measured at SPring-8 [126] and 
GRAAL [127] ; the A polarization was deduced in [127] and 
[128]. Target asymmetry T and O x and O z for jp —> K + A 
were reported in [129]. In [130], CLAS data on the spin 
transfer coefficients C x and C z were presented for both, 
■yp — > K + A and jp K + S°. 

Differential cross sections on the reaction 7p — > K°U + 
were measured by CLAS [128] , Saphir [131] and CB-ELSA/ 
TAPS [132]. For the latter data we include the determi- 
nation of the P polarization derived from an analysis of 
the S + decay. 

2.8 The reactions *yp — >■ p7r°7r° and 7p — pTv°q 

The two reactions jp — > pir ir [33,34] and jp —> pir°r) 
[35-37] are included event by event using an extended 
likelihood method. The quality of the fit can be judged 
from the description of Dalits plots. For the -fp — >• pir°ir° 
reaction we constructed Dalitz plots in m 2 a versus m 2 
with 20 x 20 bins, for the four 100 MeV 7p invariant mass 
intervals from 1350 to 1750 MeV. For the 7p — > ptjtt 
reaction, m 2 ^ versus m 2 ^ Dalitz plots were constructed 
for seven 100 MeV 7p invariant mass intervals from 1700 
to 2400 MeV. The number of bins with nonzero Monte 
Carlo events and x 2 j^bin are given in Table 2. The beam 
asymmetries [38,39] and the helicity dependence E [40] 
are included in the fit in the form of histograms. 



Both these reactions have been studied intensively, see 
[133-146] for the first and [147-149] for the latter reaction. 
For these data only selected histograms are available; they 
arc not included in our fits. 



3 Partial wave amplitudes 

A general expression for the decomposition of the two- 
particle scattering amplitude A(s, t) into partial wave am- 
plitudes A!}? (s) which describe production, propagation 
and decay of a two-particle systems with fixed total angu- 
lar momentum J, parity and (if conserved) C-parity can 
be written as: 

A( a> t)=x;^'(*)Q^U(fc)i^.::^Q^U(9) « 

where fcj are initial and g, are final particle momenta, s = 
(fcl + k 2 ) = ( Ql + q 2 ) = P 2 ,t = (h - qi) 2 = (k 2 ~ 9 2 ) 2 , 
k = (fci — k 2 )/2, q = (qi — q 2 )/2 and n = J for a boson 
system and n = J — 1/2 for a fermion one. The vertices 

Qvi...v n and Q^i„.ft n (T stands for hermitian conjugation) 
describe the transition of the system into the initial- and 
final-state particles, and depend on the total and relative 
momenta. The indices ft and ft' list quantum numbers 
of the production and decay amplitudes, e.g. isospin, spin 
and orbital angular momenta. The tensor Fj^"^ n depends 
only on the total momentum P and describes the tensor 
structure of the partial wave. It is often called projection 
operator. The formalism for construction of vertices for 
meson-baryon partial waves and projection operators is 
given in [150, 151]. For convenience we provide key formula 
for projection operators and vertices in Appendix A. 

In the case of resonance production, the total ampli- 
tude A(s, t) can be expanded into a sum of partial wave 
amplitudes multiplied by vertices, see eq. (1). Here the 
partial wave amplitudes A^f (s) provide the energy de- 
pendence of the resonance which can be parameterized, 
for example, as N/D amplitude, as K-matrix or, in the 
simplest Breit-Wigner amplitude [152]. For non- 

resonant contributions, like t and u channel exchanges, the 
situation is different. In many partial wave analyses (in- 
cluding the present one) these contributions are simply 
added to the resonant part of the total amplitude and 
the sum is used to fit the experimental data. However, 
one needs to know the contribution of t and it-exchanges 
in every partial wave if the final partial wave amplitudes 
are to be compared with results from other analyses. This 
decomposition is also required when rescattering between 
non-resonant and resonant parts of the amplitude should 
be taken into account. For the non-resonant contributions 
used in the energy dependent fits one has therefore to solve 
an inverse task: to extract partial wave amplitudes from 
the total amplitude. 

This task can be solved by using the orthogonality con- 
dition for partial wave operators. Multiplying the total 
amplitude from eq. (1) with initial and final projection 
operators and vertices and integrating over solid angle of 
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the initial and final momenta we obtain 
df2k dfi, 



= (-i) n F^:ZYl^'( s ) w n p ( k l) w n a '(9i) . (2) 

/3/3' 

where fc^ and are squared relative momenta orthogonal 
to the total momentum of the system P (see Appendix A) . 

The factor W" 13 corresponds to the on-shell one-loop 
amplitude for transition between two vertices Q^...^ n . 
can be calculated as 



It 



W r f(ki) 



TpOi l . . . OLf 



Zn = (-l) n FZl:;ZZg filUl ...g flnVn . (3) 



For meson-nucleon and jN vertices, the W"^ were calcu- 
lated in [151]. For convenience we provide the correspond- 
ing expressions in Appendix B and expressions for partial 
wave amplitudes for photoproduction of a single meson 
are given in Appendix C. 



3.1 Parameterization of the partial wave amplitudes 

In the present analysis, the partial waves at low energies 
are described in the framework of a K-matrix/P-vector 
approach. High- mass resonances (above 2.2 GeV) are de- 
scribed by relativistic multi-channel Brcit-Wigncr ampli- 
tudes. In the case of photoproduction reactions, the regge- 
ized t- and u-channel amplitudes were added to the reso- 
nant part. Then the multipolcs were calculated by solving 
eq. (2). 



3.1.1 Pion induced reactions in K-matrix approach 

The multi-channel amplitude is given by the matrix A(s) 
where the matrix element A a t(s) defines the transition 
amplitude from state 'a' to state 'b'. In eq. (1) this am- 
plitude is denoted as A^f (s) to emphasize the different 
spin-parity contributions. Now we will use the notation 
A a b{s) which identifies the initial and the final channels, 
e.g. 7 AT, 7rA, ijN, KA, irA, and omit the indices rep- 
resenting the partial wave. Scattering between different 
channels is taken into account explicitly in the K-matrix; 
the amplitude is given by 



A(s) = K (I — ipky 



(4) 



where K is the K-matrix, I is the unity matrix and p 
is a diagonal matrix of the according phase space. For 
two-particle states (for example tvN), the phase space is 
calculated simple loop diagram (see [151]). For J = 
L + 1/2, the so-called '+' states, the phase space is equal 
to 



P+( s ) 



a L 2\k\ 2L+1 k w +m N F{k 2 ) 
2L+1 V« 2m N B(L,r,k 2 



(5) 



and for '-' states with J = L — 1/2, the phase space is 
given by 



a L 2|k| 2L+1 k w +m N F(k 2 ) 
L y/s 2m N B{L,r,k 2 ) 



(6) 



where s is the total energy squared, k, is the relative mo- 
mentum between baryon and meson, k its three-vector 
component, kxo is the energy of the baryon (with mass 
mjv) calculated in the c.m.s. of the reaction. J is the to- 
tal, L the orbital angular momentum of the baryon-plus- 
meson system, and the coefficient is equal to: 



n 



2n- 1 



(7) 



The phase volume is regularized at large energies by a 
standard Blatt-Wcisskopf form B(L, r, k 2 ) with r = 0.8 fm, 
and a form- factor F(k 2 ) of the type 



F(k 2 



yl + 0.5 
A + k 2 



or F(k 2 



yl + 2.5 

A + s 



(8) 



Fits with both parameterizations yield nearly identical re- 
sults. The parameter A were taken from our previous anal- 
ysis [124,153] and fixed to 1.5 for the first parameteriza- 
tion and 3.0 for the second one. The exact formulas for 
the three-body phase volume are given in [151]. 
The K-matrix K is parameterized as follows: 



K„ 



(a) (a) 
\- 9a '9b f 

L> Ml - s + /a "' 



(9) 



where M a and <j4 are the mass and the coupling constant 
of the resonance a, and where f a b describes a direct (non- 
resonant) transition from the initial state a to the final 
state b, e.g. from irN — > AK. 

For most partial waves it is sufficient to assume that 
f a b are constants. The Sn and S31 and waves require a 
slightly more complicated structure, we use 



fab 



f (l) 
Jab 



f(2) 
Jab 



(10) 



Here the and sg b arc constants which are determined 
in the fits. In the case of the Sn wave, this more flexi- 
ble parameterization is required to describe ttN — > Ntt, 
ttN — > Nrj, and r/N — > Nrj transitions. Let us note that 
this form is similar to the one used by SAID [6] . 



3.1.2 The photoproduction amplitude 

The photoproduction amplitude can be written in the P- 
vector approach [154]. The P-vector amplitude for the ini- 
tial state 'a' photoproduction is then given by 



P b (I - lf>K)ba 



(11) 
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The production vector P is parameterized as: 

(a) (a) 

p b = T, g -^- + h (12) 



Ml 



where g^j are the photo-couplings of the resonance a and 
where non-resonant production of a final state b is de- 
scribed by contributions In general, these are functions 
of s but mostly, a constant fa is sufficient. 

The P-vector approach is based on the idea that a 
channel with a weak coupling can be omitted from the K- 
matrix. Indeed, adding to the K-matrix the ^yN channel 
would not change the properties of the amplitude. Due 
to its weak coupling, the -fN interaction can be taken 
into account only once; this is done in the form of a P- 
vector. Loops due to virtual decays of a resonance into 
Nj and back into the resonance can be neglected safely. 
A similar approach can be used to describe decay modes 
with a weak couplings. The amplitude for the transition 
into such a channel can be written as D- vector amplitude, 

A a = b a + [k{i - ipk)- 1 p] ab D b , (13) 

where the parameterization of the D-vector is similar to 
the parameterization of the P-vector: 



D, 



(a) (a) 

^ Ml - s 



(14) 



Here g^ is the coupling of a resonance to the final state 

and d b is a non-resonant production from the K-matrix- 
channel b to the final state. As in the case of the P-vector 
approach, channels with weak couplings can be taken into 
account only in their final decay, and are not taken into 
account in the rescattering. Let us note that if the final 
state is already included as one of K-matrix channels, the 
amplitude (13) reproduces the K-matrix amplitude (4). 

In cases where both, initial and final coupling con- 
stants are weak, we use an approximation which we call 
PD-vector. In this case the amplitude is given by 



A ab = G ab + Pa(I ~ ipK)' 1 pD b 



(15) 



where G a b corresponds to a tree diagram for the transition 
from state 'a' to state 



Gab = 



{a) (a) 
9 a 9 b 

Ml-s 



h a b 



(16) 



Here g\ is the production coupling of the resonance. For 
photoproduction, gi"^ = g^ holds true, and h a b is the 
direct non-resonant transition from the initial to the dif- 
ferent final channels. 



peaks originate from meson exchanges in the t-channel. 
Their contributions are parameterized as 7r, p(tu), K or 
K* exchanges. 

The most straight forward parameterization of particle 
exchange amplitudes is the exchange of Regge trajectories. 
The invariant part of the t-channel exchange amplitude 
can be written as [152] 



T(s,t) = 9l (t)g 2 (t)R{±,v,t) u=-(s-u). (17) 

Here, gi are vertex functions, and R(+, v 1 1) and R(— , v 1 1) 
are Reggeon propagators for exchanges with positive and 
negative signature. Exchanges of n and K have positive, 
p, u> and K* exchanges have negative signature. 

The p trajectory has a negative signature and the cor- 
responding propagator is equal to 



Rp{-,v,t) 



-if ot p {t) 



cos(^a p (t))r 



>(*) 



(18) 



where a p (t) = 0.50 + 0.85t. The w trajectory is identical 
to the p trajectory. The expressions for other Reggeon 
propagators used in the fit are given in Appendix D. 



4 Partial wave analysis 

4.1 Fit of the n°p and 7r+n photoproduction reactions 

The new CLAS data on the jp — > ir°p reaction are com- 
pared to our fit in Fig. 5. The % 2 contributions of this fit 
from the various channels are given in Tables 1-4. We re- 
mind the reader that we estimated additional systematic 
errors for the jp — > n°p and 7p — > ir + n differential cross 
sections; these additional errors are not shown in Fig. 5. 

Some systematic deviations between data and fit can 
be recognized in the mass region below 1800 MeV. These 
are mostly the result of discrepancies between the data. 
In Fig. 6 we show for comparison some CB-ELSA data, 
and in Fig. 7 some GRAAL data, for invariant masses 
which are close to the CLAS values. At higher energies 
the solution describes very well the new CLAS data, how- 
ever CB-ELSA data are also described with rather good 
accuracy. 

The new data on the 7^ — > 7r + n and the fit curve are 
shown in Fig. 8. Here, the total normalization factors re- 
solve rather well discrepancies at masses below 1600 MeV. 
The description of the earlier data on the jp — > n + n in 
this mass region is shown in Fig. 9. 



4.2 Photoproduction multipoles 



3.2 Reggeized meson exchange amplitudes 

At high energies, angular distributions of photo-produced 
mesons exhibit clear peaks in the forward direction. These 



We now turn to a discussion of the partial wave ampli- 
tudes. It should be stressed that the amplitudes we give 
for 77) — > pir° and 7p — > nir + are constrained by a large 
number of other reactions. This is particularly important 
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Fig. 5. CLAS data on the differential cross section for jp — > 
ir a p with current solution. Only statistical errors for the CLAS 
data are shown 
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Fig. 7. GRAAL data on the 7p — > ir p differential cross section 
at energies close to the CLAS values. Only the errors quoted 
by the GRAAL collaboration are shown. 
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Fig. 8. CLAS data on the 7p — > -K + n differential cross section 
cm with the current solution. 



Fig. 6. CB-ELSA data on the 7p — > Vp differential cross 
section for energies below 1.8 GeV. Only the errors quoted by 
the CB-ELSA collaboration are shown. 
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Fig. 9. The 7p —¥ it n differential cross section with the cur- 
rent solution. The data are taken from [75,77-82,85] 



Surprisingly, the multipoles for jp — > mr + arc in much 
better consistency. The differences in the E 2 and M 2 can 
be assigned to additional A3/2- (1940) and Z\ 3 / 2 -(2260) 
resonances introduced to fit data on ^p — > pir°r] [36,37]. 
Significantly different are the multipoles leading to 5/2~ 
states. In our fits, the E 2 and M 2 multipoles include an 
additional resonance N 5 / 2 - (2060) [153]. We note that the 
(dominant) resonance contributions are compatible with 
the Watson theorem. 

4.3 Properties of contributing resonances 



in the vicinity of thresholds. Of course, the elastic irN scat- 
tering amplitude and the pion photoproduction amplitude 
are influenced by opening new channels and the couplings 
to the new channels can be estimated from their effect on 
the scattering and photoproduction amplitudes. But this 
is rather indirect, and it is desirable to take the inelastic 
channels into account directly. 

The multipoles for 7r° photoproduction are shown in 
Fig. 12 in comparison to the SAID SP09K2700 [21] and 
MAID 2007 [22] solutions, those for jp n+n in Fig. 13. 
The errors cover a large number of fits which differ mostly 
by the parameterization of the 2-kN channel at masses 
above 1.8 GeV. For convenience of the reader, we list in 
Table 5 the lowest photoproduction multipoles and the 
corresponding partial waves. 

Most amplitudes derived within the SAID, MAID, or 
BnGa approach yield consistent results, at least qualita- 
tively. The best agreement is found for the M± amplitude 
which describes the spin flip amplitude for the photo- 
induced transition from the proton to the A resonance 
and its excitations. The A resonance is fully elastic, hence 
the agreement in the low-mass region is not unexpected. 
Even the small E^ multipoles are not inconsistent. Some 
multipoles which we discuss next show significant differ- 
ences between the different approaches. The Eq multipole 
has a similar structure in all three approaches but shows 
significant differences in detail. In the BnGa solution, the 
electric dipole transition Eq exceeds the other solutions in 
the threshold regions, likely due to a larger role of the sub- 
threshold AK + amplitude. The differences are even larger 
for the M]~ multipole; this may be not unexpected in view 
of the notorious difficulties with the l/2 + partial wave. 



Table 5. Photoproduction multipoles and partial waves. In 
general, two multipoles lead to one spin-parity wave. 



Multipoles 


Partial 


waves 


J P 


Et ~ 


S11 


S31 


1/2- 


- M- 


P11 


P31 


1/2+ 


Ef M+ 


Pm 


P33 


3/2+ 


E 2 - M 2 - 


Drs 


D33 


3/2" 


E+ M+ 


D l5 


D35 


5/2 _ 


£ 3 - M 3 " 


Fib 


F35 


5/2+ 



A large number of resonances is identified in the fits. Some 
have a strong coupling to pion photoproduction, for oth- 
ers, the product of squared photocoupling constant and 
Nir decay branching ratio is small and they contribute 
mostly to inelastic channels; their properties will be dis- 
cussed elsewhere. These latter resonances are listed in Ta- 
ble 6. They do help to improve the fit to pion photopro- 
duction but their helicity amplitudes are not well defined, 
and photo-couplings and decay branching ratios of these 
states can be varied within large limits without significant 
X 2 deterioration. All these solutions were included in the 
error estimation procedure. 

The pole position of the states, photo-couplings and 
ttN branching ratios for the states contributing strongly 
to pion photoproduction are given in Table 7. The inclu- 
sion of the new CLAS data rather notably stabilized the 
solution and improved most of the errors. The pole posi- 
tions are determined by finding zeros of the real and the 
imaginary part of the denominator in the partial wave 
amplitudes. Thus two lines are defined in the complex en- 
ergy plane. Their crossing points defines the pole position. 
The coupling constants, including the helicity amplitudes, 
are calculated as residues of the P-vector/K- matrix am- 
plitude at the pole position and are given together with 
their phases. The ttN branching ratios are calculated as 
squared residue-couplings, multiplied by the phase volume 
taken at the Breit-Wigner resonance mass. We note that 
nuclcon-meson or nuclcon-photon couplings arc defined at 
the pole position of a resonance, and are complex num- 
bers. In Table 7 we give the (complex) photon-couplings at 
the pole positions; their analogues, the helicity amplitudes 
A l j 2 and A3/2 arc defined for Breit-Wigner amplitudes, 
not for more general formalisms. The method how we de- 
rive Breit-Wigner parameters and helicity amplitudes is 
discussed below. For the Pi 3 (1720) and £> 33 (1700) reso- 
nances, the Breit-Wigner width is much larger than one 
might expect from the pole position. In the Nir channel, 
the visible width is much closer to this expectation. The 
effect is known from ao(980) which has a visible width 
of about 50 MeV in the irrj mass distribution but a much 
larger width in the KK mass distribution the width is 
much larger because of the rapidly opening KK phase 
space. In the Pi 3 (1720) and P 33 (1700) case, the phase 
space for Ntttt 3-body decays grows rapidly with increas- 
ing mass. 

Within the quoted errors, the results from the new 
solution are mostly compatible with those published in 
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W cm , MeV W cm , MeV W om , MeV W cm , MeV 

Fig. 10. The real (two left-hand columns) and imaginary (two right-hand columns) part of multipoles for the ty° photopro- 
duction. The errors are systematic and cover a large number of fits (see the text). The dashed curves correspond to the SAID 
solution SP09K2700 [21] and the dotted curves to the MAID solution 2007 [22] 




W cm , MeV W cm , MeV W cm , MeV W cm , MeV 

Fig. 11. The real (two left-hand columns) and imaginary (two right-hand columns) part of multipoles for the jp — > ir + n 
reaction. The errors are systematic and cover a large number of fits (see the text). The dashed curves correspond to the SAID 
solution SP09K2700 [21] and the dotted curves to the MAID solution 2007 [22]. 
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Fig. 12. The real (two left-hand columns) and imaginary (two right-hand columns) part of the isospin-1/2 photoproduction 
multipoles. The errors are systematic and cover a large number of fits (see the text). The dashed curves correspond to the SAID 
solution SP09K2700 [21] and the dotted curves to the MAID solution 2007 [22] 
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Fig. 13. The real (two left-hand columns) and imaginary (two right-hand columns) part of the isospin-1/2 photoproduction 
multipoles. The errors are systematic and cover a large number of fits (see the text). The dashed curves correspond to the SAID 
solution SP09K2700 [21] and the dotted curves to the MAID solution 2007 [22]. 
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Table 6. Baryon resonances included in the fit which con- 
tribute little to photoproduction of pions. 

Mass Width Mass Width 

Pn(1860) 1900±30 300±40 Pi 3 (1900) 1960±30 185±40 

Di 3 (1700) 1730±40 310±60 Z>i 3 (1875) 1870±25 150±40 

P 33 (1600) 1640±40 480±loo P 33 (1920) 1950±40 330±50 

Pis (2000) 1910±50 360±80 Z>i 5 (2070) 2065±25 340±40 

D 33 (1940) 1995±40 360±50 Di 3 (2170) 2160±35 370±50 

D 33 (2360) 2360±50 480±80 5*81(1900) 1955±30 335±40 



[33]. The largest changes to our previous solution are ob- 
served for the photo-couplings of the 5 3 i(1620) resonance 
(which was (130± 50) GeV" 1 / 2 xl0 3 in [33]), and for the 
small helicity component Ai/ 2 of the Z?i 3 (1520) (which 

was (7. Oil. 5) GeV _1 / 2 x 10 3 ). The changes are largely due 
to the inclusion of additional polarization data for 7p — > 
7r°p and jp — > ir + n. 

The new data also require Pn(1710). In [34], this reso- 
nance improved the description of the data slightly but we 
were not forced to introduce it. In the present fit, there are 
three resonances above the nucleon in the P\\ wave: the 
Roper resonance Pn(1440), the Pn(1710), and the newly 
proposed Pn(1860). 

We found two minima for the photo-couplings of the 
F35(1905) state. Both solutions are given in the Table 7 
(5 th and 6 th row, 2 nd column). The first solution corre- 
sponds well to the PDG average values, while the second 
solution has an almost vanishing hclicity-3/2 coupling. 
Both solutions reproduce single meson photoproduction 
data with the same quality, however the second solution 
provides a better likelihood for the 7p — > ir rjp reaction. 
The analysis of the high energy data on the two-pion pho- 
toproduction will help to define which solution is the phys- 
ical one. Also the forthcoming double polarization data 
will provide new and important constraints on the ampli- 
tudes. 

The pole structure of two resonances is found to be 
ambiguous. The pole of the «Sii(1650) resonance is located 
between the AK and UK thresholds. In one set of accept- 
able solutions wc found large couplings of this state to 
these channels and a complicated pole structure of two or 
even more poles close to the two thresholds. We included 
such pole positions and corresponding residues in the er- 
rors given in Table 7. The forthcoming analysis of the data 
on the itp — > KA should help to define the S\\ pole struc- 
ture more accurately, however new data on ttN — > KA and 
irN —> KS would be extremely valuable. The fit of two- 
pion photoproduction data suggests a substantial coupling 
of Pi 3 (1720) to the P>i 3 (1520)tt channel. Thus we observe 
here a double pole structure near the P>i 3 (1520)7T thresh- 
old which results in rather large errors and in difficulties to 
identify the corresponding Breit-Wigner parameters. We 
believe that the forthcoming polarization data on two-pion 
photoproduction will improve significantly the accuracy in 
the definition of Pi 3 (1720) pole structure. A more precise 
definition of the Pi 3 (1720) pole with its Nirir couplings 



Table 7. Pole position (in MeV), photo-couplings calculated 
as residues in the pole (in GeV~ 1//2 10 3 , phases in degrees) and 
corresponding Breit-Wigner parameters for states contribut- 
ing strongly to pion photoproduction (the branching ratios are 
in percents). The PDG values are given in parentheses. For 



P 3 5(1905) two solutions are given. 


State 


Sn(1535) 


Sii(1650) 


Re(pole) 
-2Im(pole) 

A 1/2 (w) 

Mbw 

Pbw 

Pn-JV/P 


1510±25 (1510±20) 
140±30 (170±80) 
90±25 /0°±45° 

1535±20 (1535±10) 
170±35 (150±25) 
35 ± 15 (45 ± 10) 


1670±35 (1655±15) 
170±40 (165±15) 
65±30 /28° ± 15° 

1680±40 (1658±12) 
170±45 (165±20) 
50 ± 25 (78 ± 18) 


State 


Pi i (1440) 


S 3 i(1620) 


Re(pole) 
-2Im(pole) 

A 1/2 (-f P ) 
Mbw 

rBw 

Pirjv/P 


1370±4 (1365±15) 
193±7 ( 190±30) 
-48±12 /-58°±20° 
1440±12 (1445±25) 
335±50(325±125) 
60 ± 6 (65 ± 15) 


1596±7 (1600±10) 
130±10 ( 118±3) 
62±10 /-0°±20° 

1625±10 (1630±30) 
148±15 (143±8) 
23 ± 5 (25 ± 5) 


State 


Pn(1710) 


P 33 (1232) 


Re(pole) 
-2Im(pole) 

A 1/2 {~ip) 
A 3/2 (~/ P ) 
Mbw 
Few 

PjtJv/P 


1708±18 (1720±50) 
200±20 (230±150) 
24±8 /-20°±60° 

1725±25 (1710±30) 
200±35 (150±100) 
12 ±6 (15 ±5) 


1211±1 (1210±1) 
100±2 ( 100±2) 

-136±5/-17°±5° 
-267±8 /-3°±3° 

1230±2 (1232±1) 
112±4 (118±2) 
100 (100) 


State 


Pi 3 (1720) 


P 33 (1700) 


Re(pole) 
-2Im(pole) 

A 1/2 (~f P ) 
A 3/2 (~f P ) 
Mbw 
J^bw 

PrAr/P 


1660±35 (1675±15) 
360±80 (190±85) 
140±50 /-35°±25° 
110±50 /10°±35° 
1770±100 (1725±25) 
650±120 (225±75) 
14 ±5 (15 ±5) 


1650±30 (1650±30) 
275±35 (200±40) 
160±45 /35°±12° 
165±40 /40°±18° 
1780±40 (1710±40) 
580±120 (300±100) 
16 ±7 (15 ±5) 


State 


Pi 3 (1520) 


P 35 (1905) (sol.l) 


Re(pole) 
-2Im(pole) 

A 1/2 (-yp) 

A^ 2 ( JP ) 

Mbw 

Fbw 

P.iv/P 


1512±3 (1510±5) 
110±6 (112±7) 
-30±6 /15°±10° 

130±6/6°±5° 
1524±4 (1520±5) 
117±6 (112±13) 

57 ± 5 (60 ± 5) 


1800±15 (1830±5) 
300±20 (282±18) 
28±10 /-35°±15° 
-42±12 /-25°±15° 
1890±25 (1890±25) 
335±30 (335±65) 
12 ± 3 (12 ± 3) 


State 


Di S (1675) 


P 35 (1905) (sol.2) 


Re(pole) 

-2Im(pole) 

A 1/2 (j P ) 

^ 3/2 (7P) 

Mbw 

Fbw 

PttJv/P 


1650±5 (1660±5) 
143±7 ( 138±12) 
20±4 /-6°±6° 
24±8 /-6°±6° 
1678±5 (1675±5) 
177±15 (148±18) 
37 ± 5 (40 ± 5) 


1805±15 (1830±5) 
310±20 (282±18) 
47±10 /-30°±12° 
0±3 

1850±20 (1890±25) 
345±30 (335±65) 
12 ± 3 (12 ± 3) 


State 


Pis (1680) 


P 37 (1950) 


Re(pole) 
-2Im(pole) 

A 1/2 (~f P ) 
A^ 2 ( JP ) 
Mbw 
Fbw 

PjtJv/P 


1672±4 (1673±8) 
114±12 (133±12) 
-12±6 /-45°±30° 

130±8 /0°±10° 
1685±5 (1685±5) 
117±12 (130±10) 

66 ± 8 (68 ± 3) 


1882±8 (1880±10) 
262±12 (240±20) 
-81±8 /-15°±12° 
-93±8 /-15°±15° 
1928±8 (1933±18) 
290±14(285±50) 
44 ± 8 (40 ± 5) 
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Table 8. Helicity amplitudes A 1/2 and A 3/2 for iV* and A* from this work, from SAID08 [21], from MAID07 [22], from the 
Giefien model [23,24] and estimates from Ref. [1]. 



Resonance 




A 1/2 (GeV 


"^xlO 


) 






A 3/2 (GeV" 1/2 xl0 3 ) 






BnGa09 


FA08 MAID07 


GieBen 


PDG 


BnGa09 


FA08 


MAID07 


Giefien 


PDG 


Sn(1535) 


90±15 


100.9±3.0 


66 


95 


90±30 












Si i (1650) 


60±20 


9.0±9.1 


33 


57 


53±16 












Pi i (1440) 


-52±10 


-56.4±1.7 


-61 


-84 


-65±4 












Pn(1710) 


25±10 






-50 


9±22 












Pi3(1720) 


130±50 


90.5±3.3 


73 


-65 


18±30 


100±50 


-36.0±3.9 


-11 


35 


— 19±20 


£>i 3 (1520) 


-32±6 


-26±1.5 


-27 


-15 


-24±9 


138±8 


141.2±1.7 


161 


146 


166±5 


£>is(1675) 


21±4 


14.9±2.1 


15 


9 


19±8 


24±8 


18.4±2.1 


22 


21 


15±9 


Pis(1680) 


-12±6 


-17.6±1.5 


-25 


3 


-15±6 


136±12 


134.2±1.6 


134 


116 


133±12 


S 3 i(1620) 


63±12 


47.2±2.3 


66 


-50 


27±11 












P 33 (1232) 


-136±5 


-139.6±1.8 


-140 


-128 


-135±6 


-267±8 


-258.9±2.3 


-265 


-247 


-250±8 


P>33(1700) 


160 ±45 


118.3±3.3 


226 


96 


104±15 


160±40 


110.0±3.5 


210 


154 


85±22 


P 35 (1905) 


28±12 


11.4±8.0 


18 




26±11 


-42±15 


-51.0±8.0 


-28 




-45±20 


or: 


(48±12) 










(0±3) 










P 37 (1950) 


-83 ±8 


-71.5±1.8 


-94 




-76±12 


-92±8 


-96±8 


-121 




-97±10 



may also help to resolve the discrepancies in the determi- 
nation of its A 3 / 2 helicity amplitude when different anal- 
yses are compared. 

The Review of Particle Properties lists Breit-Wigncr 
parameters and real helicity amplitudes. In the K-matrix 
approach, a photo-produced resonance is described by a 
P-vector (eq. 12). Even if the photo-coupling constant 

9^n ^ s rGai i t ne photo-coupling at the resonance position, 
calculated as residuum of the amplitude Pb at the pole, 
will in general be complex. To allow for a comparison 
with other determinations, we define helicity amplitudes 
by the following procedure: a Breit-Wigncr amplitude is 
constructed with an adjustable mass and a width which is 
parameterized as a sum of all partial widths, Pigf- The 
total widths is scaled with one parameter. This scaling 
parameter as well as the Breit-Wigner mass are adjusted 
to reproduce the pole position of the P-vector amplitude. 
The Breit-Wigner helicity amplitudes A x / 2 and A%/ 2 are 
defined by the condition that the residues of the Breit- 
Wigncr photoproduction amplitude reproduce the mag- 
nitude of the original residues of the P-vector/K- matrix 
amplitude. 

In Table 8 we compare our results on A x / 2 and A 3 / 2 f° r 
N* and A* with previous determinations of these quan- 
tities. These real helicity amplitudes are given in Table 8 
and compared to values obtained by SAID, MAID, and 
the Giefien model, and to the values listed by the PDG 
[1]. First we notice that our errors are much larger than 
those given by FA08, MAID and Giefien do not give any 
errors. We believe that the FA08 systematic errors are 
underestimated: the impact of variations in the couplings 
to inelastic channels can hardly be tested using only re- 
actions with Ntt in the final state. The errors we quote 
are not statistical errors; those are small. Our errors are 



derived from a large number of fits changing the number 
of resonances, switching on and off couplings to inelastic 
channels, using different start values for the fits. 

For most resonances, reasonable consistency between 
the different analyses is found. In particular the helicity 
amplitudes for photoproduction of the Roper resonance 
from SAID, MAID, and BnGa are fully consistent (the 
Giefien result is a bit higher) even though mass, width, 
and Ntt decay branching fractions differ somewhat. BnGa 
and SAID, e.g., find, respectively, 

M (MeV) r (MeV) JWAot 

BnGa 1440 ± 12 335 ± 50 0.60 ±0.06 
FA08 1485 284 0.79 

In our analysis, the Roper resonance is fully constrained: 
from three of the four reactions, irN elastic scattering, 
7P — ¥ Nir, ir~p — > p7r°7r°, and 7p — > p7r°7r°, the ampli- 
tude for the forth reaction can be predicted. Hence we are 
particularly confident that these results are correct. 

We comment briefly on further differences. The PDG 
result for the A 1/2 amplitude of (53±16) GeV" 1 / 2 x 10 3 
for producing Sn(1650) was driven by the 1995 VPI re- 
sult (69±5)GeV" 1 / 2 xl0 3 [155] and by the small value 
(22±7)GeV~ 1/2 xl0 3 obtained in [47]. The most recent 
FA08 analysis gives (9.0±9.1) GeV" 1 / 2 x 10 3 , a value which 
is much smaller and which is not confirmed here; we find 
(60±20) GeV~ 1/2 x 10 3 , in close agreement with the Giefien 
result. Part of the discrepancy with FA08 is certainly due 
to the Sn(1650) branching ratio to the irN channel; in 
FA08 this is fixed to be 100% while we find (50±25)%. Of 
course, photoproduction defines only the product of the 
helicity and irN couplings. 

Possibly related are the differences in the helicity am- 
plitudes for Pi3(1720). Our value for A 1 / 2 is compati- 
ble with the new FA08 analysis and in conflict with the 
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value quoted by the PDG. Incompatible with all other 
determinations - even in the sign - is our value for the 
helicity amplitude for Pi 3 (1720) production. Also 
the Giefien results are at variance with the other deter- 
minations. Clearly, more data are required to resolve this 
discrepancy; the results from double polarization experi- 
ments carried out at present in different laboratories will 
very likely be decisive. 

There is the possibility, that the discrepancies in the 
properties of Pi3(1720) have a physical origin. In [156] 
it was found that data on the reaction ep — > e'pn + Tr~ 
could be described only when resonance parameters were 
drastically changed with respect to published results, or 
when a new resonance in the Pi 3 wave was introduced. 
Apparently, the Pi3(1720) properties are different in Nir 
and in Ntttt; this might be a hint for the presence of a 
close-by state in the same partial wave. 



5 Summary 

We have presented results from a partial wave analysis 
on a large variety of different reactions, from ttN elastic 
scattering to photoproduction of multibody final states. 
The main emphasis of this paper was devoted to a deter- 
mination of the electric and magnetic multipolcs leading 
to the production of neutral or charged pions in photo- 
induced reactions off protons. The multipoles are mostly 
consistent with previous analyses but a few significant 
discrepancies call for clarifications. The analysis provides 
masses, widths, and helicity amplitudes for several known 
resonances. Masses and widths and the irN partial decay 
widths of all resonances agree very well with established 
values. Only the photocoupling of the Pi3(1720) resonance 
differs remarkably from PDG and from the values found 
in a recent analysis of the CLAS collaboration. This dis- 
crepancy may be a further hint for the conjecture [156] 
that the Pi3(1720) resonance may have a more compli- 
cated structure than usually assumed. 
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qi + q2 and two relative momenta fc^~ and qj^ which are 
orthogonal to the total momentum: 



:{qi - Q2)u9L, 



9 [ii/ 



9uu 



PyPu 
p2 



(19) 



The tensor Fj£-"£ n depends only on the total momen- 
tum P (s = P 2 ) and describes the tensor structure of 
the partial wave. It can be calculated as a product of two 
polarization tensors „ summed over possible polar- 
izations: 



(20) 



For every set of indices, P^ 1 "^™ satisfies the properties of 
the polarization tensor: it is symmetrical over permutation 
of any two indices, traceless and for n > is orthogonal to 
the total momentum of the system. It usually normalized 
by the condition: 



(21) 



and is often called projection operator: its convolution 
with another tensor by one set of indices results in a tensor 
which obeys the symmetry properties of the corresponding 
partial wave. 

In the case of a fermionic system, P^ 1 "^" can be writ- 
ten in the form 

pPi-M~ _ ( 1 \n Vs + P (~,fti— H„rp£i~4n (~)0i.-.p n (oo\ 

Here, (y/s + P) corresponds to the numerator of a propa- 
gator describing a particle with J = 1/2 and n = J— 1/2 
(s/s = M for the stable particle). We define 



m£i..4„ n+l, n 
1 



(24) 



We introduced the factor l/(2y/s) in the propagator which 
removes the divergency of this function at large energies. 
For the stable particle it means that bispinors are normal- 
ized as follows: 



u(k N )u(k N ) = l , 



u(kisr)u{kN) = — — — .(25) 



polarizations 



2m 



Appendix A: The structure of the fermion prop- 
agator 

We consider scattering of two particles with momenta k\ 
and in the initial and q\ and q2 in the final state. There 
are three independent momenta. It is convenient to choose 
the total four-momentum of the system P = k± + k2 = 



Here and below, k = 

The boson projection operator Of\ 
ing properties: 



;f " has the follow- 



1 y.i KJ vi_...v n — Jr Vj^'v 1 ...v„ 



Of} 



— n ■■■**» — f) 



nCl-Cr. f\a 1 ...oi n 
^ai ...a n v\ ...v n 



(-l) n O%;;;£ 



(26) 
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For the lowest states, 

= 1, OS=g^, 



I / 2 \ 

^i^2 2 [9fi 1 u 1 9fj. 2 u2^ 9f J . 1 u 2 9fj. 2 u 1 ^ "^9ii 1 fj. 2 9v 1 iy 2 

.(27) 

For higher states, the operator can be calculated using the 
recurrent expression: 

1 / " 

^vi...v n n 2 \ ' yywj KJ v i _...vj-iVj+\...v n 



(2n- l)(2rc-3) 



E„-L n ^ f}/ti...ftj_i/t i+ i.../^_ I /^ + i...Jt„ \ /r>o\ 
y j Ui/lji'!/fcI/ m , ^!/l...I/Jt_l!/|i + l...!/ m _ll/ OT +l...!/ n I ' 



»<3 
fc<m 



The tensor Fj£"£ n has all orthogonality properties of the 
tensor 0£ 1 ""£ n plus orthogonality to the 7-matrix: 



IPi-T Vt...V n r u 1 ,,,U„ IV} — U 



(29) 



The pscudoscalar meson-nucleon vertices for the par- 
tial wave with spin J have the form: 



CU„ = *757^+!k (?X<7i) , (30) 



where n = J— 1/2 and is the bispinor of the baryon. 

The '+' and '-' indices describe two sets of the partial 
waves with relation between orbital momentum L and the 
total spin J as J = L+l/2 ('+' partial waves) and J = L — 
1/2 ('-' partial waves). The '+' set of vertices describes the 

partial waves with J p = i , | + , | ... and the second 

S6t J — 2" 5 2 '2 ■ • ■ • 

In the case of virtual photons there are, for every par- 
tial wave with J > 1/2, three j*N vertices; for real pho- 
tons, only two of them are independent [151]. In the LS 
formalism these vertices correspond to spin 1/2 and 3/2 
of the photon-nucleon system. For '+' states the vertices 
arc (following the ordering in [151]): 



Q { £l« n = "(fei)7,ns4S ) 1 ...a„_ 1 (fc ± )^ n ^ . (3i) 

where u(k±) is bispinor of the initial nucleon and is the 
photon polarization vector. 
For '-' states we have: 



)(!- 



[if 1 



(32) 



The orbital angular momentum operators for L < 3 



<*>=z^<x> 



Fig. 14. Diagram representation of eq. (2) 



Y(2) - 3 ( jU-L jU-L _ 1 t.2 ± 



M1M2 / ' 



A"(3) = £ 



jU-L jU-L jU-L 



The operator X^^ 1 .^ can be written as a series of 
products of metric tensors and relative momentum vec- 
tors. The first term is proportional to the production of 
relative momentum vectors k't, other terms correspond 
to the substitution of two vectors by a metric tensor with 
corresponding indices: 



n+l 



k 2 

k X k X k 1 - k 1 - k X - -± 

^Si>a>3-">» 2n+l 



e 9^ n k H + e su k » n k L + ■ 

j 4 / n 

I )(2n— 1) \ ^ 9vnt9ii jt i m Y\. k f*i 



(2n+l)(2n-l) 



i<k,j<m 



(34) 



Appendix B: Contribution of the loop diagrams 

For the irN vertices we have [151]: 

, m N + fcio 



= (-1) 



n "n ll.|2n' 



2n + 1 277i 7v 



n + l 



2uin 



(35) 



and for the jN vertices (in the case of photoproduc- 
tion): 



Wi 11+) = (-l)"2^Hk| 2 " mjV + fcl ° , 
2n + 1 2mjq 

w ^ . ( _ ir _^_(ii±l) |kr ii| 



+ fc 



in 



AT 



^i 13+) = (-^"^rlkl 5 
2n + 1 

for the '+' states and 



, m N +k w 
2mjv 



(36) 



X (0) 



W (ll-) = / 1 J a "+l | k |2n+2 TO Af + fclO 

n+l 2mj\r 
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ttH33-) ( n n a n -i(n + l) , 2ra _ 2 mjy + fcio 
1 j (2n+l)(2n-l)' 1 2tojv 



13- 



<-i)-£l|k 



2n TO JV + fc 10 



(37) «„■(+,!/,*) 



for the '-' states. The 7iV vertices in this representation 
are not orthogonal to each other, and to extract partial 
waves one needs to solve a 2 x 2 system of linear equations. 



Appendix C: Single meson photoproduction am- 
plitude 

The general structure of the single-meson photoproduc- 
tion amplitude in c.m.s. of the reaction is given by 

t~ t~ i _ (crk) _,_ (era) 

k q k q p q 2 



(38) 



where q is the momentum of the nucleoli in the -kN chan- 
nel and k the momentum of the nucleon in the 7-ZV channel 
calculated in the c.m.s. of the reaction. The <r, are Pauli 
matrices and oj, u>' are non relativistic spinors of initial 
and final states correspondingly 

If T% are known, e.g. from the t or u channel exchange 
amplitudes calculated in the c.m.s. of the reaction, the 
partial wave amplitudes can be obtained as 



A (i±) _ f dz {i±) 
~ I 2 m ' 



(39) 



where z is the cosine of the angle between initial and final 
relative momenta and vectors D^- 1 ^ are equal to 



D< 1+ ) = 

D (2+) = 

DC 1 -) = 

D (2-) = _ 



1 1 








l-z 2 


(0,0, ■ 






n - 


Hi / 



n+l , 
P' 



(1-^+1 

Vi+l)(n+2) 

n P' 



Kn+l&n+l 
l-Z 2 
Kn-ia„-l|k| : 



(n+l)' (n+l)(n+2) y 
p p (l- 2 )^ + i 

r n-\-l i r ri i 



(n + l)(n+2),0 



o o ni " +1 P' 

' ' (n+2)' " 



(40) 



Here P n = P n (z) are Legendre polynomials and P' n = 
dP n (z)/ dz. 

Using the multipole decomposition of the A^' ampli- 
tudes given in [150] one can obtain the standard expression 
for the projection of the total amplitude into multipoles. 



Appendix D: Reggeon propagator parametriza- 
tion 

In this section we give the expressions for Reggeon prop- 
agators used in the fit. 



The propagator for pion exchange has the form 

-if a„(t) / ,, \ 



e 2' 



sin(fc(t))r(^ + l) V^o 



(41) 



where a v (t) = —0.014 + 0.72t is a function defining the 
trajectory, Uq is a normalization factor (which can be taken 
to be 1). The _T- function is introduced in the denominator 
to eliminate the additional poles at t < 0. The propagator 
for Kaon exchange is given by 



R K (+,v,t) = 



sm(^a K (t))r ( 



2 



CMC (t) 



,(42) 



where a K (t) = -0.25 + 0.85i. 

The propagator for K* exchange is identical to the p 
exchange propagator but has «k« (t) = 0.32 + 0.85t. 
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